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A model is described that relates the temperature coefficient of the optical absorption spectrum 
of pure water to the frequency derivative of that spectrum and two parameters that quantify the 
dependence of a peak's amplitude and its position on temperature. When applied to experimental 
temperature coefficients, it provides a better understanding of the process than the analysis currently 
in use. 
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I. INTRODUCTION 

Understanding the temperature dependence of optical 
absorption takes on increasing importance as more ac- 
curate spectra become available. A high quality absorp- 
tion spectrum from the infrared to the blue for water 
at 22° C can be obtained by combining the data of Kou, 
Labrie, and Chylek[l[ with that of Pope and Fry Kou 
et al. studied the range from 666 to 2400 nm with a 
Fourier-transform infrared spectrometer. They derived 
absorption coefficients from the ratio of transmittances 
obtained with different path lengths. This exactly can- 
cels reflection losses at the air- window and window- water 
interfaces and absorption by the cell windows. Only 
transmittance values between 20% and 60% were used 
to avoid saturation and base line uncertainties. Absorp- 
tion coefficients from 380 to 727.5 nm were obtained by 
Pope and Fry using a novel integrating cavity absorption 
meter (IC AM). In the ICAM, optical power loss from a 
water-filled highly-reflective cavity is measured as a func- 
tion of wavelength. Scattering, which is important at the 
lower end of this wavelength range, does not cause a loss 
of power, hence the ICAM measures absorption. 

The absorption coefficient of water decreases by some 
six decades from the infrared through the blue. The in- 
frared spectra in Fig. 1 from Kou et al. |l[ show a general 
exponential decrease vs frequency that is modulated by 
various peaks, valleys, and shoulders indicating absorp- 
tion by discrete states or groups of states. The amplitude 
of these features decreases relative to the background 
continuum with increasing frequency. The use of fre- 
quency (in wavenumbers) emphasizes the approximately 
uniform spacing of these features. 

The vibrational level spectrum of an isolated water 
molecule can be described in terms of harmonics or com- 
binations of three normal modes: a symmetric and an 
asymmetric stretch, and a bending or scissors mode. 
Single levels are broadened by rotational structure. In 
condensed phases, formation and breaking of hydrogen 
bonds further broadens levels and may introduce addi- 
tional modes. The result is not a continuum for low exci- 
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FIG. 1: Absorption spectra of water, supercooled water, and 
ice from Kou, Labrie, and Chyleky|. Oscillator assignments 
are shown in parentheses above the peaks, see text. The ratio 
of the 22° spectrum to the -8° spectrum is shown in the lower 
part of the figure. Vertical dashed lines mark inflection points 
on the high frequency side of each 22° peak. 



tations: Spectra in Fig. 1 retain various discrete features. 
Patel and Tarn Q have identified absorption peaks in the 
liquid with the band structure of two anharmonic oscilla- 
tors. The first was based on a single (degenerate) stretch- 
ing mode. The second, a parallel series of states, was up- 
shifted by a bending excitation. We will refer to these as 
the stretching band and the combination band, respec- 
tively. The frequency in wavenumbers of levels in terms 
of stretching and bending numbers, (n s , rib), is given by: 



u(n s ,0) — n s (3620 — 63n s ] 
v(n s ,l) = v(n a ,Q) + 1645. 



(1) 
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The peaks and shoulders in Fig. 1 are labeled with their 
(n s ,rib) assignments. 

The 22°C absorption spectrum of Pope and Fry [2[ at 
has a minimum in the vicinity of 420 nm. Recent data 
from Abe and the Super-Kamiokande Collaboration^ 
at 13°C, and from Griskevich, Renshaw, and SmyQ at 
15°C, suggest that this minimum is appreciably deeper 
and shifted to a higher frequency. The present paper is 
an outgrowth of a study of scattering and absorption by 
water in the vicinity of the minimum. When compar- 
ing various measurements, it was desirable to separate 
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the influence of temperature from those due to other sys- 
tematic effects. The classical treatment of temperature 
dependence is to fit the absorption spectrum to an appro- 
priate set of Gaussians and then to determine a unique 
temperature coefficient for each peak from experimental 
data. An example of this sort of analysis due to Jonasz 
and Fournier [6| is described in Sec. 3. Measurements 
of Buiteveld, Hakvoort, and Donze [7|, of Pegau, Gray, 
and Zaneveld Q , and of Langford, McKinley, and Quick- 
enden Q indicate that absorption peaks in the visible 
not only shift upward, but also move to higher frequen- 
cies with increasing temperature. The model developed 
below predicts that shifts of peak position may result in 
apparent positive and negative temperature coefficients 
for a given peak. A variety of experimental observations 
are a natural consequence of this model. 



II. THE MODEL 

Ratios of absorption coefficients at 22°C to those at 
— 8°C are plotted in the lower part of Fig. 1. Although the 
individual absorption spectra show considerable struc- 
ture, the ratio takes a simple form; a series of cusps where 
the ratio is greater than one, separated by regions where 
it drops below unity. Vertical dashed lines in Fig. 1 indi- 
cate that the cusps correspond with inflection points on 
the high frequency sides of peaks in the 22° C spectrum, 
not their centers. This suggested a connection between 
temperature dependence and the shape of the absorption 
spectrum. 

Consider part of an absorption spectrum, a(v,T), in 
the vicinity of some spectral feature located at frequency 

where the absorption coefficient is clq . We assume that 
a is separable: 

a(y,T) = a o (T)f(v-v (T)). (2) 

Differentiating Eq. 2 with respect to T and substituting 
dfJdvQ = —df/dv leads to: 

1 da 1 da Q 1 da dv Q 
a dT ao dT a dv dT 

The temperature coefficient 1 is the sum of two terms. 
The first is the trivial "classical" interpretation: If a tem- 
perature change of +1°C increases the intensity of an ab- 
sorbing peak by +1%, then the temperature coefficient 
is +1% per °C across the whole peak. 

The second term on the right side of Eq. 3 has more 
interesting properties. If the region of interest contains 
a peak, da/dv is zero at the peak position. The temper- 
ature coefficient will change sign in scanning across that 



1 We will use "temperature coefficient" for the fractional tempera- 
ture coefficient, reserving "absolute temperature coefficient" for 
da/dT. 



peak. That maximum temperature changes match the 
inflection points on the high frequency side of the peaks 
in Fig. 1, indicates that dvo/dT is positive in the infrared. 
Experimental data presented in Sec. 3 indicate this is true 
in the visible as well. As a consequence, temperature co- 
efficients will be negative on the low frequency side of a 
peak and positive on the high side. 

If the region of interest is a decreasing continuum, the 
derivative is negative throughout, and temperature co- 
efficients will be positive. Consider then the progressive 
transition from a spectral region dominated by individ- 
ual peaks to one of a decreasing continuum. In the first 
region, temperature coefficients should show a bipolar os- 
cillation, with some upwards bias due to the underlying 
continuum and any contribution from the classical term 
of Eq. 3. Moving up in frequency, the relative contribu- 
tion of peaks decreases and that of the continuum in- 
creases. The amplitude of the oscillations will decrease, 
and the upward bias will increase. At some point, the 
negative lobes will be balanced out by the bias, and an 
absorption peak will appear to have only a positive tem- 
perature coefficient peak. However, this will not be cen- 
tered at the peak position but on the inflection point on 
its high frequency side. With further increases in fre- 
quency, the envelopes of the upper and lower coefficient 
oscillations will converge to the continuum value. 



III. EXPERIMENTAL DATA 

The logarithmic frequency derivative, (l/a)(da/du), 
is compared with the fractional temperature coefficient, 
{X/a){da/dT), in Fig. 2. The derivative of the 22°C spec- 
trum in Fig. 1 is shown as a hollow curve in Fig. 2(a). 
The solid curve is from the absorption spectrum of Pope 
and Fry [2|, also at 22° C. Peak positions of the stretching 
and combination band absorptions are shown by (n s , n.&) 
labels. Each of the oscillator states is marked by a sharp 
negative derivative spike on its high frequency side. Pos- 
itive excursions are less pronounced due to, in part, the 
overall decreasing trend of the spectrum. These deriva- 
tive oscillations damp down in moving into the visible. 
Negative excursions corresponding to each state persist, 
but the positive excursions transform into a small nega- 
tive background between the larger negative excursions. 
Pronounced peaks and valleys have given way to an expo- 
nentially decreasing absorption spectrum where specific 
absorbing states appear only as shoulders. Note that the 
amplitude of negative derivative peaks associated with 
(n s , 1) transitions decreases more rapidly than those for 
(n s ,0). They are decidedly less pronounced in the visi- 
ble. This suggests that the dependence on frequency is 
different for the two bands. Positions of the derivative 
minima are projected on to Fig. 2(b) by vertical dashed 
lines. 

The dependence of the temperature coefficient on fre- 
quency is shown in Fig. 2(b). The hollow curve was ob- 
tained from the data of Kou et a/.[l|. The solid curve 
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FIG. 2: Comparison of the frequency derivative of the water 
absorption spectrum at 22° C with the temperature coefficient 
of that spectrum. Hollow curves in (a) and (b) are from Kou 
et al. The solid curve in (a) is from Pope and FryQ]. The 
solid curve in (b) is based on the measurements of Langford 
et al. fj|. A dashed curve connects it to the hollow curve. Cir- 
cles show experimental values from Pegau et aZ.[j|, squares 
those from Trabjerg and Hojerslev The heavy dotted 

curve is from an analysis by Jonasz and Fournier [6j. Thin 
dotted curves show the approximate envelopes of the temper- 
ature coefficient oscillations. Positions of frequency derivative 
minima are shown by vertical thin dashed lines. 

is based on the extensive measurements of Langford, 
McKinley, and Quickenden These authors studied 
the temperature dependence of decadic absorptivity, A, 
over the wavelength range from 550 to 900 nm and for 
temperatures between 15 and 60° C. At any wavelength, 
A was observed to vary linearly with T . Their published 
decadic dA/dT values were converted to da/dT, then to 
(l/a)(da/dT) using the absorption coefficients of Kou et 
al. and of Pope and Fry. A spline interpolation joining 
the two measurements is shown as a heavy dashed line 
in Fig. 2(b). 

Buiteveld et a/.@] have derived da/dT values for wave- 
lengths from 400 to 800 nm (12,500 to 25,000 cm" 1 ) from 
absorption spectra obtained between 2.5 and 45°C. They 
noted that da/dT was nearly constant at 0.0012 m -1 per 
°C below 580 nm. Pegau et al. have concluded there 
was a systematic bias of 0.0011 in that value. We sub- 



tracted that bias prior to converting the Buiteveld et al. 
da/dT values to (1/ a) (da/dT). These were in agreement 
those of Langford et al. for wavelengths longer than 580 
nm. They have not been included in Fig. 2(b) because of 
uncertainties due to the baseline correction, 

Absolute temperature coefficients from Pegau et aZ.[|[ 
were converted to the (1/ a) (da/dT) values which are 
shown as circles in Fig. 2(b). With the exception of the 
point at 850 nm (11,760 cm -1 ), there is good agreement 
with the curve based on Kou et al. and Langford et al.. 
Four points due to Trabjerg and Hojerslev [lfj, squares 
in Fig. 2(b), are also in agreement with that curve. 

Jonasz and Fournier 6] have reported a detailed anal- 
ysis of the absorption spectrum of water and its tem- 
perature coefficient in the range from 380 to 800 nm 
The absorption spectrum could not be reproduced with a 
ten-member basis-set of Gaussians corresponding to the 
stretching and combination band states (4,0) to (8,1). 
Five additional Gaussians were needed to obtain "impres- 
sive" agreement with the experimental spectrum. Jonasz 
and Fournier then used available data to derive frac- 
tional temperature coefficients for each of the 15 Gaus- 
sians. This classical approach is equivalent to setting 
dv$/dT = in Eq. 3. Their derived spectrum of tem- 
perature coefficients, the dotted curve in Fig. 2(b), is a 
reasonable fit to the experimental data for frequencies be- 
low 15,500 cm -1 where the extra Gaussians were added. 
In particular, it reproduces the negative values to the left 
of the (4, 0) peak and the near zero values in the vicinity 
of the (4, 1) position. Agreement is poorer at higher fre- 
quencies. Oscillations of their curve in that region are a 
consequence of their ascribing temperature coefficients of 
0.45% per °C to the (n s ,0) absorptions and 0.20% per °C 
to the (n s ,l) absorptions in their original ten-member 
basis set. The oscillations would have continued down to 
12,500 cm -1 (800 nm) were it not for the five extra Gaus- 
sians. One positioned near the (4, 0) peak had a larger 
positive coefficient. The other four located between the 
(4, 0) and (5, 0) peaks had large negative values. The 
physical significance of unusual properties of these phan- 
tom states is not clear. 

Several features of Fig. 2 support the derivative 
model. The one-to-one correspondence between fre- 
quency derivative minima and temperature coefficient 
maxima is striking and implies that dv§ / dT in Eq. 3 is 
positive. The nearly symmetric, two-lobed form of the 
temperature coefficient spectrum in the infrared provides 
strong evidence that frequency shift is a major contribu- 
tor to temperature dependence. The region of near zero 
temperature coefficients at ~ 15,000 cm -1 is interpreted 
as being due to a balance between the negative lobe of 
the (5,0) absorption and and the positive contribution 
from the decreasing continuum. The (4, 1) absorption 
makes a small contribution in this valley. 

Approximate upper and lower envelopes of the tem- 
perature coefficient oscillations for the (n s , 0) absorptions 
are shown as thin dashed curves in Fig. 2(b). The damp- 
ing of oscillations with increasing frequency, and the the 
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approach to a limit of ~ +0.3% per °C which is shown as 
a heavy dashed line, are consistent with the change from 
discrete absorption peaks in the infrared to a decreas- 
ing continuum in the visible. Note that the 0.3% per °C 
applies to absorption only. Because scattering has little 
if any temperature dependence, the fractional tempera- 
ture coefficient for attenuation is reduced by the ratio of 
absorption to attenuation. For example, scattering and 
absorption coefficients are equal at 425 nm Q and the 
fractional temperature coefficient for attenuation is then 
0.f5%pcr°C. 

IV. CONCLUSIONS 

The temperature dependence of optical absorption by 
water has a complex dependence on frequency related 
to properties of the absorbing states or continuum. The 
present model, which connects the temperature depen- 
dence to the frequency derivative of the absorption spec- 



trum, provides a framework for understanding this de- 
pendence. The composite curve in Fig. 2(b) based on 
Kou et al. and Langford et al. with the dashed interpo- 
lation and extrapolation, provides practical temperature 
corrections for absorption coefficients. The model sug- 
gests that analyses of the type which fit an absorption 
spectrum with Gaussians should use the derivative of the 
Gaussians to fit the temperature coefficient spectrum. A 
continuum background must be included in both. 
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